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In high density QCD the hadron production stems from decay of mini-jets that have the trans- 
verse momenta of the order of the saturation scale. It is shown in this paper that this idea is able 
to describe in a unique fashion both the inclusive hadron production for ^/s > 546 GeV including 
the first data from LHC and the deep inelastic scattering at HERA. Recently reported data from 
ALICE, CMS and ATLAS including inclusive charged- hadron transverse-momentum and multiplic- 
ity distribution in pp collisions are well described in our approach. We provide predictions for the 

J5 ' upcoming LHC measurements. 

^' 

G^ ' I. INTRODUCTION 

(N' 

The first LHC data [l|-3 on inclusive hadron production call for a theoretical understanding of these processes 

,^ ■ based on QCD. At first sight the inclusive hadron production is a typical process that occurs at long distances where 

Qh' one has to use the non-perturbative methods of QCD. Therefore, the field of long distance processes seems to be a 

f-^ ■ relevant subject to the domain of high-energy phenomenology with the main ingredients: soft Pomeron and secondary 

Q ' Reggeons. Such phenomenology is able to describe inclusive hadron production data (see Ref. Q and references 

^S^ , therein) but cannot be considered satisfactory since both soft Pomerons and Reggeons cannot be explained in terms 

of QCD ingredients; quarks and gluons. It should be also mentioned that the increase with energy of the average 

transverse momentum of the produced hadron observed experimentally [3, |3| cannot be explained in the Reggeon 

J> ' approach. 

T—i However, high density QCD j6l-[l^ leads to a completely different picture of inclusive hadron production. In this 

CO approach the system of parton (gluons) at high energy forms a new state of matter: Color Glass Condensate (CGC). 

^5 ■ In the CGC picture, at high energy the density of partons pp with the typical transverse momenta less than Qs reaches 

a maximum value, pp ex l/a^ 3> 1 (as is the strong coupling constant). Qs is the new momentum scale (saturation 

I/"} ■ momentum) that increases with energy. At high energies/small Bjorken-x, Qs ^ p where p is the scale of soft 

^D , interaction. Therefore, as (Qs) ^ 1 and this fact allows us to treat this system on solid theoretical basis. On the 

other hand, even though the strong coupling as becomes small due to the high density of partons, saturation effects, 

the fields interact strongly because of the classical coherence. This leads to a new regime of QCD with non-linear 

features which cannot be investigated in a more traditional perturbative approach. 

In the framework of the CGC approach the secondary hadrons are originated from the decay of gluon mini-jets with 
the transverse momentum equal to the saturation scale Qs{x). The first stage of this process is under theoretical control 
C^ and determines the main characteristics of the hadron production, especially as far as energy, rapidity and transverse 
momentum dependence are concerned. The jet decay, unfortunately, could be treated mostly phenomenologically. 
However, we can hope that the phenomenological uncertainties would be reduced to several constants whose values 
will be extracted from the experiment. 

Actually, such a description has passed the first check with the experimental data: the KLN paper [ij] explains the 
main features of inclusive hadron production in heavy ion-ion and hadron-ion as well as proton-proton collisions |14l | 
at RHIC. In this paper we wish to improve the KLN approach by introducing two new elements: the probability to 
find gluon with fixed transverse momentum that describes the deep inelastic scattering (DIS) data and that satisfies 
the Balitsky-Kovchegov [9|, lllj non-linear equation; and a different description of inclusive hadron production at low 
transverse momenta of gluons. Over all success of our description indicates universality of the saturation physics 
which can be further tested at LHC and future collider experiment. 

In the next section we discuss the fc^-factorization and main formulas that we use. In particular, we consider the 
interrelation between the color dipolc scattering amplitude and the unintegrated gluon density that follows from the 
recent development of high density QCD |15| . An important improvement here to the previous works based on the 
KLN approach is the explicit inclusion of the impact-parameter dependence of the saturation scale. Section HI is 
devoted to comparison with the experimental data and to discussion of various predictions for higher LHC energies. 
As a conclusion, in Sec. IV we highlight our main results and predictions for LHC. 
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II. INCLUSIVE GLUON PRODUCTION IN HIGH DENSITY QCD 

The gluon jet production in hadron-hadron collisions can be described by /ct-factorization given by [li 

da 



where a:i^2 = [pt I \fs)&^^ : and pt and y is the transverse momentum and rapidity of the produced gluon jet. 0^* are 
the probability to find a gluon that carries Xi fraction of energy with fcr transverse momentum and Cp = {N"^ ~ l)/2A'c 
is the SU{Nc) Casimir operator in the fundamental representation with the number of colors equals Nc- 

For a proof of /cj-factorization see Ref. [15l| and also Refs. [16h2(j| which confirm the former proof^. We need to 
recall that the proof for the fct-factorization was given for the scattering of a diluted system of partons, say for virtual 
photon, with a dense one. Our main idea is that we have gluon saturation for proton-proton scattering or in other 
words, we are dealing with interactions of two dense systems of partons (gluons). Therefore, the fct-factorization has 
to be considered here as an assumption. It should be noticed that the proof given in Refs. [15l - l20| shows that the 
fct-factorization is valid in the situation where two scales of hardness: the transverse momentum of the produced 
gluon (pt) and the saturation scale are both larger than the scale of the soft interaction (/i). For dense-dense system 
scattering we have actually three scales: pr and two saturation scales. However, only for the kinematic region where 
both xi and X2 are small and for px which is smaller than both saturation scales we have to make an assumption 
about fc(-factorization. In other cases that one of the saturation scales is small, we are dealing with diluted-dense 
system scattering. We believe that the fcr-factorization is currently the best tools at our disposal for the processes 
considered in this paper. 

The unintegrated gluon density (I)q I xi; fc^ ] and color dipole-proton forward scattering amplitude N {xi, tt', b) are 
related in a very specific way [15[. This relation reads as follows 
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a;.; kr) = —j77^ / dHd^rre^'^^-''^ V| 7V^' {y, = ln(l/xO; rr; b) , (2) 
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(2^)3 
with 

N^^{y, = \n{l/x,);rT;b) = 2 N {y, ^ ln{l/x^);rT;b) - N^ {y, ^ \n{l/x^);rT;b) , (3) 

where Nq' [yi = ln(l/xi); r^; fo) is the dipole-hadron [hi) forward scattering amplitude which satisfies the Balitsky- 
Kovchegov equation. In the above, tt denotes the transverse dipole size and b is the impact parameter of the 
scattering. 

Eq. ([3|) looks very natural at large Nc- Indeed, for the color dipole amplitude in the Glauber form N = 1— exp (— f7/2) 
(51 is the opacity), equation Eq. ^ leads to Nq = 1 — exp {~Vt) as it should be for the scattering of the two dipoles 
of the same sizes. We recall that a colorless gluon-probc just creates such two quark-antiquark dipoles, and the Nq 
is directly related to the gluon density. 

Substituting Eq. ^ in Eq. ([T]), and after analytically performing some integrals, we obtain [1^ 



da 2Cf 1 



f dHd'Bd^fre:^^-^'^ VlN^' (j/i = ln(l/a;i);rT; 6) V^, iV^^ [y^ = ln(l/a;2); tt; \b - S|) 



dyd^PT as(27r)4 p^ 

(4) 

In the above equation, B is the impact parameter between center of two hadrons and b is the impact parameter of 
the produced mini-jet from the center of the hadron, see Fig. [TJ 

A. Choice of color dipole scattering amplitude 



As it can be seen from Eq (I2I4|) . we need here an impact-parameter dependent color-dipole forward amplitude. 
We will show later that the inclusion of the impact-parameter is very important in our approach and should not be 



^ Ref. :2l] states that Eq. JT} is not correct. Unfortunately, there is no discussions in the paper why their result is so different from the 
other pubhshed papers. However, Braun has recently shown that Ref. I21I actually leads to the fct-factorization |22|| . 




FIG. 1: Mini-jet production in hadron-hadron collisions in the transverse plane. The impact-parameter between two hadrons 
is B. 

ignored. The dipole-proton forward scattering amplitude N {Y;r;b) (with Y = ln(l/a;)) can be in principle found 
by solving the perturbative nonlinear small- a; Balitsky-Kovchegov (BK) [3, |Tl| or Jalilian-Marian-Iancu-McLerran- 
Weigert-Leonidov-Kovner (JIMWLK) [l^] quantum evolution equations. Unfortunately, numerical solution to these 
non- linear equations in the presence of the impact-parameter is very challenging [23| and is not yet available. Moreover, 
a numerical solution does not give us the full control on the phenomenological parameters that have been used and 
we certainly lose the transparency and simplicity of physical interpretation if we rely only on the numerical solutions. 
Therefore, we choose a different approach to the solution of the BK equation that was suggested in Ref. [2J]. First, 
we recall that the BK equation predicts the geometric scaling behavior [25|, namely the amplitude N {Y; r; b) is not 
a function of three variables but it is a function of only one variable Z^ = r^ Ql{x; b) ( N {Y; r;b) = F (Z)) where 
Qs{x]b) is the saturation momentum^. We also know Q the behavior of the scattering amplitude deeply in the 
saturation region (Z ^ 1) 

N{Y;r;b) = 1 - cxp f- -^^2!^ In^ zV (5) 

V 2(1 - 7cJ / 

where ^(7) is the BFKL kernel 

a; (7) = asxh) = «.{2V(1) - ^(7)) - VXI-7)}, (6) 

with a notation dig = asNc/i^. In above, we define ^^{x) — dh\V{x)/ dx and r(a;) is the Euler function. The parameter 
7cr is the solution to the following equation 

dx ilcr) ^ X ilcr) ,^. 

d^cT 1 - 7cr' 

In Ref. [26| a solution was found for the entire kinematic region for a simplified BFKL kernel, namely, instead of 
Eq. ^ , the following kernel was used, 

r i for Z = rQsix) < 1 ; 

^(7) = «s< (8) 

[ ^ for Z = rQ,{x) > 1; 

which describes only leading twist contribution to the full BFKL kernel of Eq. ^ . The lesson from this solution is 
very instructive: for r^Ql{x) < 1, the amplitude N satisfies the DGLAP (BFKL) linear evolution equation with the 



^ Notice that here we assumed that the geometric scaling is also vaUd in the presence of impact-parameter dependence of the saturation 
scale. It should be stressed that the proof of the geometric scaling behavior 25] could be easily generalized to the case of the scattering 
amplitude that depends on the impact-parameter b. In the analytical solution of Ref. J2m] which gives the theoretical basis for the 
chosen parametrization of the dipole-amplitude here, the 6-dependence is taken into account, therefore, this solution gives a theoretical 
example of the general proof. 



boundary condition N{Y; r; b) ^ Nq = Constant for r^ = 1/Q'^{x) while for r^Ql{x) > 1 we have solution that has 
the form of Eq. ([5]). Using these general features of the solution we choose the model suggested in Ref. [23| which 
improves the earlier studies on this line [2J, [2^ . In this model the color dipole-proton forward scattering amplitude 
is given by 

r ^^(|^2(,. + ^in(i)) for2 = rg.(x) <2; 
N{Y;r;b) = I (9) 

[ 1 - cxp{-A\n^{BZ)) for Z = rQs{x) > 2; 

where the saturation scale Qsix] b) (denoted by Qsix) for brevity) is given by 

Q.(.;., .(^)*»p{-j^^^JL_}^ (10, 

As we have already mentioned Eq. (J9l) as well as Eq. (fTO|) has the form of the solution to the BK equation at a 
fixed QCD coupling. For Z < 1 the effective anomalous dimension 7s + —^ In (-§:) with 7s = 1 — jcr follows from 
the BFKL (and DGLAP) equation in the vicinity of the saturation line (see Ref. [^ for the detailed derivation). 

For the leading order BFKL kernel with frozen QCD coupling the parameters of Eq. ^ and Eq. ^TU\\ have the 
following values 

l-j„ = 0.63; A = Qs ^ " = 4.88as; k = ^" y'^''' ^ 9.9. (11) 

1 - Icr X'hcr) 

The parameters A and B can be found from matching of N and its logarithmic derivatives at Z = 2 while A'o and 
BcGC remain fitting parameters. 

Generally speaking, the model given by Eqs. (I9ll0p can be viewed as an approximation to the solution of the BK 



equation. However, because the 6-dependent numerical solution to the BK equation is not yet available j23[, we are 
doomed to resort to such an approximation. This model differs from other saturation models on the market since it 
apparently incorporates all known properties of the exact solution to the BK equation including the 6-dependence of 
the scattering amplitude (see Ref. [26,]). 

The advantage of Eqs. (|9ll0p is that these equations give the possibility to take into account the next-to-leading 
order (NLO) corrections. Two features of the non-linear low- a; equations can be calculated in the next-to-leading 
order using the kernel of the linear equation: the energy behavior of the saturation scale [6i, [32, |33[ and the behavior 
of the solution deeply in the saturation domain 26]. It has been shown that the NLO correction to the BFKL 
equation (and therefore BK equation) are large and it changes considerably the value of A from A « 0.9 to A sa 0.3 
for ckg = 0.2 29, 3 G|. The value of 7s in Eq. ([SJ is also affected by the NLO corrections as well as by the running 
QCD coupling [29l - [3l| . It is therefore generally believed that the higher order corrections to the NLO BK equation 
should be important. The actual calculation of higher-order corrections to these non-linear evolution equations still 
remains as a challenge. Since the general behavior of the amplitude Eq. (9) will remain unchanged after inclusion 
of higher-order corrections, we effectively incorporate the higher-order corrections by taking the value of parameters 
A,7s, A'o and Bcgc obtained from a fit to the DIS data at low Bjorken-a; x < 0.01 [23] ■ Therefore, the saturation 
model that we use here gives also a good description of the HERA data at low-z. In order to simulate the behavior 
of gluon density at large x — )■ 1, we product the unintegrated gluon density with (1 — x)^ as prescribed by quark 
counting rules jSj]. This factor stems from the correct description of the HERA data on DIS. 

B. Physical observables 

The rapidity distributions of the mini-jets can be calculated using Eq. ([1]), 

'-^^=^ - m^ fd^PT j^ [E,. ©] , (12) 

dri cJnsd J dy d^pr 

where 77 is the pseudorapidity and h[r]] is the Jacobian which takes account of the difference between rapidity y and 
the measured pseudo-rapidity 77 ^13i] , 

M.,PT)=^^2i^_, (13) 

^^^ +sinh^ 



where rrijet is the mass of mini-jet. One also has to express rapidity y in Eq. ([T]) in terms of pseudo-rapidity 77. This 
relation is given by 



-, I w — ^ + smh 77 + sm?] 

y{v,PT) = 7:lnj ^, — }. (14) 
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The distribution Eq. ([T]) refers to the radiated gluons with zero mass while what is actually measured experimentally 
is the distribution of final hadrons. We therefore should make an assumption about hadronization of gluons which 
is entirely non-perturbative process that has to be modeled in any approach due to lack of understanding of the 
confinement of quarks and gluon in QCD. However, it is well-known that the general assumption about hadronization 
leads to the appearance of mass of the mini-jet which is approximately equal to rri^^f ~ 2/ipT (see Ref. [131] ) where 
/i is the scale of soft interaction. The mini-jet mass rrijet effectively incorporates the non-perturbative soft pre- 
hadronization in the pseudo-rapidity space. Accordingly, one should also correct the kinematics every where in 



Eq. ([T} due to the presence of a non-zero mini-jet mass, namely replacing pT — > \/Pt + ^et i^ ^i' -^2 and also in the 

denominator of l/px- One cjin see that Eq. ([1]) has infrared divergence at pt —^ for the kinematic region kr ^ pt 
when rrijet = 0. In Ref. [131 it was suggested to integrate over fc^ < Pt- The reason is that such an integration 
reproduces the factorization formula at large pr ^ 1^ for the DGLAP evolution. However, as we explained above it 



is more natural to replace pT by ^^/px + ^et ^^ ^'i- © which consequently also regulates the denominator due to 

the presence of a non-zero mini-jet mass (the appearance of such mass is the general property of the hadronization 
processes) . 

In Eq. ([T2|) we do not take into account the fragmentation of the produced gluon (mini-jet) into hadrons. We rely 
on the principle of Local Parton-Hadron Duality (LPHD) [35, 36:] namely the form of the rapidity distribution will not 
be distorted by the jet decay and only a numerical factor will differ the mini-jet spectrum from the hadron one. We 
believe that it is better to use the LPHD scheme than to deal with the fragmentation's functions for which we have 
no theoretical justifications at low pt- It should be stressed that the same idea has been used in the KLN approach 
which describes the rapidity distribution of heavy-ion collisions data in a wide range of energies. This idea has also 
worked perfectly in e+e~ annihilation into hadrons |35l.l36l]. 

We should stress that the value of inelastic non-singlet diffractive (NSD) cross-section Oj^d cannot be calculated in 
our approach and has to be taken from the soft interaction models such as in Refs. [37|, [S^]. The NSD cross-section 
ansd is defined as ansd = o'tot — o'ei — o'sd — <7dd where aei, <Jsd and Udd are the cross sections of elastic, single and double 
diffraction respectively. However, the experimental data on add is very limited [33], cr^d is measured with rather large 
errors [40, |4l[ and even for the total cross-section utot |4l| we have two values at the Tevatron energies [42| . Therefore, 
we should stress that in this way we can only predict da/dy rather than dNch/dy. In order to overcome this problem, 
here we choose a different strategy: the physical meaning of ansd in Eq. (J12p is the area of interaction which can be 
calculated in our approach. Indeed, using Eq. (j4]) one can calculate the average impact parameter for the inclusive 
production of the mini-jet 

b'.et) = (15) 

^d'^bd^Bd^rriy + 16- S]^) e*'^""^""' V^N^^iyi = ln(l/xi); tt; 6) V|, 7V^' (y2 = Hl/x2)]rT;\h-B\ 
j£^J d%d^Bd\Te^^T--T v^TV^i (yi =ln(l/a;i);7-T;&) V| A^^^ (772 - ln(l/x2); ry; Jfc - B]) 

The NSD cross-section a^sd is then equal to the average interaction area upto a constant (TjysD — -^^^(^?et)- ^^^ 
pre-factor M will be determined and discussed later. We should draw the reader attention that such a picture for the 
inelastic cross-section corresponds, in a sense, to the geometric-scaling behavior of the scattering amplitude. Indeed, 
the high-density QCD deals with the partonic wave-function of a fast hadron which describes a coherent system of 
partons (quarks and gluon). At high energy the coherence of partons is destroyed during a short time, and the partons, 
distributed as in the wave function, are produced. These partons contribute to the inelastic cross-section. The elastic 
(diffractive ) cross-section corresponds to a rare event where the target does not destroy (or destroyed only partially) 
the coherence of the gluons in the wave- function (see for example Ref. [43]). The geometric-scaling behavior as well 
as the saturation phenomenon, in general, means that partons are distributed uniformly in the transverse plane in the 
wave-function of a fast hadron in a such way that the wave-function generates a uniform distribution of the produced 
partons after the interaction with the target. Therefore, the NSD (inelastic) cross-section is proportional to the area 
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occupied by partons. Actually, such a view on the inelastic cross-section was suggested in the KLN approach [13[ but 
for nucleus-nucleus and hadron-nucleus collision. Therefore, we generalize this approach to hadron-hadron scattering. 
We believe that if the LHC data at higher energy will support this idea, it will be a strong argument in favor of 
the saturation approach. The relation (tnsd = crtot ~ cTei — (^sd — crdd shows the obvious fact that the prediction for 
elastic and diffractive scattering are much more complicated and less transparent in the saturation approach. This is 
well-known fact at least for diffractive production .44|. 

The average transverse momentum of the mini-jet is defined in the usual way: 

(Pjet,T) = j dvh[v]j d^PT\pT\^^^[Eq.m]/ j dilh[if^j d?pT^^^^[Eq.m- (16) 

The advantage of this quantity is that it can be calculated without usual uncertainties associated with the soft 
interaction physics. The average transverse momentum of the jet can be directly related to the saturation scale via 
Eqs. (I1I9I16P and it has the following simple form at large Qs ^ rUjet, 

where the parameter Q is of order of one and takes into account the contribution of integrals in Eq. (|16p for px > Qs- 
In order to calculate the transverse momentum of hadrons which is measured experimentally, we need to recall that 
^hadron.T = ^Pjct.T + Pintrinsic.T ^^^ich leads to 



'^hadron.T/ y '-^^jet.T/ + '^intrinsic.T/ > \^°) 

where z is the fraction of energy of the jet carried by the hadron. (^intrinsic t) ^^ ^^^ average intrinsic transverse 
momentum of the hadron in the mini-jet. In other words, this is the transverse momentum of the hadron in the 
mini-jet that has only longitudinal momentum. 

In the framework of the LHPD, the pr spectrum of the produced hadron is equal to 

^%dron ^ [drjhl^j]^ —^ \Eq. Q with pj^t ^ = pr/z] , (19) 

d^PT J CTnsd d-qd^p-„i^rp L i""^'^ J 



jet,T 



where in the above px is the transverse momentum of the produced hadron. 

In the CGC scenario, the gluon saturation scale is proportional to the density of partons (see Refs. [13, [Tj])- The 
parton density is proportional to the multiplicity and, therefore, we can use the following expression for the saturation 
momentum in the event with the multiplicity of the hadrons n: 

Ti 

Qs(x)^Qs{n-x) = --Q,{x), (20) 

{n) 

where (n) is the average multiplicity that has been measured in inclusive production without any selection related to 
multiplicity and Qs{x) is the saturation scale for inclusive hadron production or Qs{n =< n >;x). Using Eq. p7)) 
again one can relate the saturation scale at a given multiplicity to the average transverse momentum of the produced 
mini-jets at large Qs ^ rUjet, 



III. COMPARISON WITH THE EXPERIMENTAL DATA AND PREDICTION FOR HIGHER 

ENERGIES 

In the derivation of the fct-factorization it was assumed that the strong coupling as is a constant. As a generalization, 
in Eq. ([I} we replace as by asipr), where pr is the transverse momentum of the mini-jet and in Eq. ^ we also 
replace 1/as by l/as{Qsixi)) where Qsixi) is the saturation scale in hadron hi. This seems to be the most natural 
way of introducing the running coupling which still preserves the form of Eq. Q apart from the over-all factor outside 
of integrals which now depends on kinematics. Indeed the inclusion of running strong-coupling leads to improvement 
of our description. For the running strong coupling as , we employ the same scheme as used by the KLN approach 
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FIG. 2: Right: shows the average impact parameter of the produced mini-jet ^^^^^ given by Eq. (|15p as function of energy. Left: 
The comparison with the experimental data and prediction for dNch/dy using Eq. (|12p with anad ~ MTv{b'j^t) for |?7| < 2.4. The 
curves are normalized by data at ^/s = 546 GeV, see the text for the details. The experimental data are from Refs. [ ll. I2I. l4q . 
The error bars on the UA5 and ALICE data points are statistical. We show only systematic errors for the CMS data points. 



[13|, namely we use the leading-order running coupling with smooth freezing below the virtuality Q^ « 0.8 GeV at 
the value of af^ « 0.5. This is in accordance with many evidence from jet physics which indicates that the QCD 
coupling may stay reasonably small, al^ = 0.4 -^ 0.6 in the infrared region 45]. 

The impact-parameter dependence in our formulation emerges from the employed impact-parameter dependent 
saturation scale, see Eqs. (|9|10p . In this model, the profile of the saturation scale in the proton is assumed to be a 
Gaussian. It is difficult to interpret the parameter Bcgc in Eq. (jlOl) in terms of proton size due to the dipole size 
r and rapidity Y dependence of the anomalous dimension. Nevertheless, in order to have a intuitive picture, one 
may take 2Bcgc as relative average squared transverse radius of the proton. The value of Bcgc = 7.5 GeV~^ was 
obtained as a fit in order to describe the slope of t-distribution of diffractive processes at HERA [23| , which in turn fix 
the normalization of the color dipole-proton cross-section. In Fig. [5] (right), we show the average impact parameter 

of jet (6^gj) from center of the hadrons. Notice that for obtaining (??e(), the over-all coefficient in Eq. P^ will 



The 



b%, 



is about 2. 5 Bcgc and it slightly increases with 

2fipT 



13| since the typical transverse 
The saturation scale in 



be dropped out and we are left with no free parameter, 
energy. 

The mass of mini-jet rrijet is proportional to the sat 
momentum of the mini-jets is the saturation scale Qs and fj, is the scale of soft interaction, 
the CGC-b model Eq. ([9]) changes slowly with energy. For our interested range of energy considered in this paper at 
midrapidity 77 = and px = I GeV for the central collisions 6 = 0, we have Qs ~ 0.6 -^ 0.8 GeV. Taking the scale of 
soft interaction equal to pion mass fi « m.^ = 0.14 GeV, we have rujet ~ 0.4 -^ 0.5 GeV. We will first assume a fixed 
value for the mini-jet mass rrijet = 0.4 GeV. To estimate the effect of the mini-jet mass, we will later consider a case 
with a different value for rrijet- 

In order to obtain the multiplicity distribution of hadrons in pp collisions from the corresponding mini-jets pro- 
duction cross-section Eqs. (|lll2p we have to fix some unknown parameters. First, based on the gluon-hadron duality, 
the rapidity distribution of hadron and radiated mini-jets can be different by a factor C. Second, although the kt- 
factorization incorporates the small-x evolution taking into account the higher-order gluon scatterings and non-linear 
gluon recombination effects, nevertheless given that we resort to a phenomenological color-dipole model, there might 
be still some extra contributions which are missed in our formulation. The discrepancy between the exact calculation 
and our formulation can be then effectively taken into account with a extra K-ia,ctoT. Finally, in order to obtain 
the charged-particle multiplicity, we should divide the mini-jet cross-section with non-singlet diffractive cross-section 

which as we already discussed is obtained via cTnsd = Mirl 5^g( ) with a new unknown dimensionless parameter M. For- 
tunately, these three unknown pre-factors C, K and M appear as a product and can be reduced to only one unknown 
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FIG. 3: Right; Energy dependence of the charged hadrons multiplicity in the central region of rapidity 77 = in pp collisions. 
The theoretical curve (Saturation model LR) is our prediction coming from the saturation model for the NSD interactions. 
The band indicates about 2% theoretical error. The total theoretical uncertainties is less 6% at high energies (see the text for 
the details). We also show the KLN prediction [Uj] with the same error band as ours. Left: Our prediction for the energy 
dependence of the average transverse momentum of charged hadrons. The CMS data ^ points and the theoretical curves 
in the left panel are for |r;| < 2.4. The experimental data are from Refs. [2. H l46r SQj . The experimental error bars indicate 
systematic uncertainties. 



parameter which will be determined with a fit to the experimental data for the charged particle multiplicity 
midrapidity for the lowest energy considered here ^/s = 546 GeV. Therefore, we obtain -^ = 2.32 at ^/s ■ 



— i — at 

ari 

546 GeV. 

We assume that this over- all normalization factor is energy-independent. We expect that the energy-dependence of 
the normalization factor to be proportional to 1 -I- 0{\/ln{l/x)). Then for higher energy ^/s > 546 GeV, we do not 
have any free parameters in our calculation and our results may be considered as predictions of the model. Notice 
that we have employed a color-dipole model that its free parameters was obtained from a fit to the HERA data for 
xb < 0.01 and Q^ G [0.25,45], therefore our formulation is less reliable at lower energies (now used here). In Fig. [5] 
(left), we show the charged multiplicity distribution for pp collisions at various energies. Our model gives a good 
description of all available data for y/s > 546 GeV including the recently released data from ALICE [l| , CMS [21 and 
ATLAS [i| at 0.9 and 2.36 TeV. We also show our predictions for the LHC energies at 7, 10 and 14 TeV. It is seen that 
as the energy increases the peak of rapidity distribution at forward (backward) becomes more pronounced. This effect 
has been also observed in Ref. [5l| where it was shown that the rapidity dependence of the invariant cross-section for 
both identified hadrons and direct photon has a peak at forward rapidities and this peak will be further enhanced by 
saturation effects [51| . 

In Fig. 13] (right) we show the charged-hadron pseudorapidity density in the central region 77 = as a function of 
center-of-mass energy in pp collisions. Notice that since our prescription is valid only for the NSD interactions we 
do not show the corresponding data for the inelastic event selection. We have also shown recently reported charged- 
particle pseudorapidity density from ALICE 3 at 7 TeV in \ri\ < 1 for inelastic collisions with at least one charged 
particle in that region (denoted by INEL > 0). Again this point is out of the scope of our calculation and we did not 
expect to explain it. 

The main source of possible theoretical error in our calculation are due to the uncertainties associated with assuming 
a fixed value for the mini-jet mass for all energies and the uncertainty in value of energy-independent normalization 
factor KC/M obtained from a fit. The value of mini-jet mass is controlled by the saturation scale and as we already 
discussed, it can be rrijet < 0.65 GeV for our interested range of energy here. Notice that the saturation scale in our 
model varies very slowly with energy. The upper limit of the theoretical uncertainty band in Fig. [3] (right) corresponds 
to a higher mini-jet mass rrijet — 0.5 GeV. The experimental systematic and statistical errors in the data point taken 
for fixing the normalization also induce uncertainty in the value of pre- factor KC/M obtained from a fit. This error 
is included in the band shown in Fig. [3] (right) and is less than the uncertainties coming from modeling the mini-jet 
mass. Over all we expect less than 6% theoretical error in our calculation at higher energies. 
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FIG. 4: The differential yield of charged hadrons for \rj\ < 2.4. The experimental data are from CMS [2| at 2.36 TeV for 
Irj] < 2.4. We show also our theoretical predictions for 7 and 14 TeV with < z >— 0.5 and rujet ~ 0.4 GeV. The experimental 
error bars shown are systematic and statistical errors added linearly. 



Our approach improves saturation based (KLN approach) calculation [l3| in several ways, including: we used a 
correct relation between the unintegrated gluon-density and the forward dipole-nucleon amplitude Eqs ()2|3p in the 
fct-factorization Eq ([T]). As it is seen this relation is not a simple Fourier transform of the dipole-amplitude which is 
commonly used in literature and also depend on the impact-parameter. The impact-parameter dependence in these 
equations is not trivial and in principle should not be assumed as an over-all factor. We then employed an impact- 
parameter dependent saturation model which was obtained from a fit to low Bjorken-x HERA data. In this sense, 
we had no freedom in modeling the saturation physics compared to the KLN approach. Moreover, since we have an 
impact-parameter formulation here, we could calculate the average relative interaction area at higher energies and 
thereby could also determine the relative increase of the NSD cross-section. It should be recalled that in the KLN 
approach the information about ansd was taken from the models for the soft high-energy interactions which is alien 
to the saturation approach. In both approaches, lower energy data for pp was used to fix the overall normalization 
factor. Therefore, we expect that the discrepancies between our predictions and the KLN to be more pronounced at 
higher energies. This is indeed the case as it can be seen in Fig. [3l it is seen that the KLN prediction underestimates 
the multiplicity at higher energies. 

The average transverse momentum of charge hadrons can be obtained from Eq. ([T5)) . In Eq. p^ . the average 
intrinsic transverse momentum of hadron has a purely non-perturbative origin and is due to the finite-size effect of 
hadrons. We take (^intrinsic t) squal to the pion mass, the scale of soft-interaction /i = to^ throughout this paper. 
In order to obtain the average transverse momentum of charge hadrons, we need also to know the value of the average 
momentum fraction of mini-jets carried by the hadrons (z). It is seen from Fig. [3] (left) that an average value of 
(z) = 0.48 -^ 0.5 is remarkably able to describe the average transverse momentum of charge hadrons in a wide range 
of energies. Our theoretical curves and CMS data are for the range \ri\ < 2.4. One may also estimate the value 
of (z) from the fragmentation functions, having in mind that the (z) for mini-jets in parton-hadron duality picture 
is not necessarily the same as the corresponding average of fragmentation momentum of the produced gluons in the 
parton model. Nevertheless, employing recently developed AKK08 fragmentation functions [53 for charged hadrons 
production from a gluon, one obtains (z) = 0.5 on average over low pr within the range of 1 < j>t [GeV] < 2 (AKK's 
fragmentation is valid only for Q > 1 GeV). In order to further test the validity of the value (z) « 0.5 for the mini-jets, 
we show in Figs. HI [5jtop panel) our predictions obtained from Eq. ([T9|) for the differential yield of charged hadrons 
in the range jryj < 2.4 and at various |?7| bins for ^/s = 2.36 TeV. The experimental data are recently reported from 
CMS collaboration 2] . It is seen that our results is in quite good agreement with experimental data. We recall again 
that the pre-factor in Eq. ()19|) is the same as what we already fixed with experimental multiplicity data at low-energy 
a/s = 546 GeV at 77 = 0. Therefore, we have no free parameters in obtaining the theoretical curves in Figs. I4I5I (top). 
In Figs. I4I5I (top), we have also shown our predictions for ^/s — 7 and 14 TeV. The fact that our model reasonably 
works at low pr (for y/s = 2.36 TeV) is due to the fact that the saturation scale is rather large at low px, for px ~ iti-w 
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FIG. 5: Upper panel: The differential yield of charged hadrons in various |ryj bins for ^/s — 2.36 TeV. The experimental data 
are from CMS ^. We also show our predictions for 7 TeV and 14 TeV with < z >= 0.5 and nijet ~ 0.4 GeV . The experimental 
error bars shown are systematic and statistical errors added linearly. Lower panel: The differential yield of charged hadrons for 
Irj] — 0.1 for two different value of mini-jet masses nijet- The inserted plot in the lower panel figure shows the charged hadrons 
multiplicity again for two values of mj^t for the same energy. 



we have Qs ~ 1 GeV in the central rapidity region. Notice that the LPHD in the simplified form that has been used 
here, is less reliable at higher pT and one should then somehow model the fragmentation of mini-jets into hadron. 

In Fig. [SI it is seen a peculiar peak of the charged hadrons production rate at low pr ~ 0.2 GeV. Actually the 
appearance of such a peak is expected in our formulation. Notice that from Eq. (J12p the differential yield of charged 
hadrons has a form '^ ^ 



drjdpi 



■2 ,,^^2 '2 ■^{xi,X2,Pt) where J^ is an analytic function. At pr — trivially we have 
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FIG. 6: The average transverse momentum of charged hadrons as a function of the number of charged particles for events with 
rich > 1 within the kinematic range pr > 500 MeV. The experimental data are from ATLAS for y's = 0.9 TeV and \r\\ < 2.5 
[3|. The theoretical curves was obtained for \r]\ = and with the same kinematic constraint pr > 500 MeV at various energies 
for two value of (z) = 0.48, 0.5 corresponding to the dashed and the solid lines, respectively. We only show the systematic 
experimental errors. 



d^JV 



(z) the spectra is a monotonically increasing function of px and for px > rrijet {z) it is 

i-jetiz) ~ 0.2 GeV since 



Tind^ - 0' for PT < mjet\ _ 

decreasing due to the denominator. The position of the peak is then approximately at px — "ttlj^ 

we have (z) = 0.5 and nijet = 0.4 GeV. This simple picture is consistent with the CMS experimental data [3| shown 

in Fig. [5] (top). 

In order to see more clearly the effect of the mini-jet mass rTijet , in Fig- E] (down) we compare the differential yield of 
charged hadrons calculated with two different values for the mini-jet mass rrijet — 0.4 and 0.8 GeV. We also show the 
multiplicity distribution in the inserted panel in Fig. [5] As we already pointed out, the mass of mini-jet is controlled 
by the saturation scale. Obviously from the saturation scale in our model, rrijet = 0.8 GeV is too large. Therefore, it 
is not surprising that the description of experimental data for both multiplicity and spectra worsened for such a large 
mini-jet mass. Nevertheless, it is obvious from Fig. [S]that the position of the peak moves to a higher pT for a larger 
mini-jet mass. Note that the CMS experimental data i^ at y^ = 2.36 TeV for the average transverse momentum of 
charged hadrons can be reproduced with (z) = 0.37 when rrijet = 0.8 GeV. Again the position of the peak in spectra 
is consistent with simple formula pt — mjet{z) ~ 0.3 in accordance with the full calculation shown in Fig. [5] Notice 
that in our model calculation shown in Fig. [5] (top), the position of the peak persists at various rapidities bin (and 
energies) since we have taken a fixed rrijet for simplification. To conclude, a precise measurement of the differential 
yield of charged hadrons at low px for higher energies at LHC will provide valuable information about the mini-jet 
mass and its connection with the gluon saturation. 

In Fig. [S] (top) , we also showed our theoretical predictions for 7 and 14 TeV with a fixed (z) = 0.5 and rrijet = 0.4 
GeV. As we already explained due to the possible increase of mini-jet mass at higher energies, the position of peak 
may slightly move to higher pr within 0.2 < pT[GeV] < 0.3 at -y/s — 14 TeV. 

In Fig. [51 we show the average transverse momentum of charged hadrons as a function of the number of charged 
particles for events within the kinematic range px > 500 MeV. The experimental data are from ATLAS for y/s — 0.9 
TeV [3| . The saturation scale at various multiplicity is given by Eq. (|20p where (n) can be conceived as a normalization 
and its value is taken to be the charged multiplicity at midrapidity rj = for a given center-of-mass energy (shown in 
Fig. [3] (right)). In order to implement in our calculation the experimental kinematic constrain px > 500 MeV on the 
measured events, we impose that (^intrinsic t) ^ '^^^ MeV. The (^intrinsic t) ^^^ ^ purely non-perturbative origin 
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and can be of order of hadron mass. To this end, we take (Pjjitrinsic t) ~ ("^p ~^ 'Tife)/2 where the mass of p and k 
mesons are nip — 775 MeV and mk — 497 MeV, respectively. In Fig. [HI we show (pr) for two values of (z). It is seen 
that our model is able to give a very good description of the ATLAS data. We have also shown in the same plot, our 
predictions for the higher LHC energies. 

The general behavior of the theoretical curves shown in Fig. |3] (left) and Fig.lHlfor the average transverse momentum 
of the produced hadrons is in accordance with a simple formulas given in Eqs. (|17l2ip showing a clear connection 
between the gluon saturation and the measured transverse momentum of charged hadrons. 

IV. CONCLUSION 

In high density QCD the main source of hadron production is the decay of gluon mini-jets with the transverse 
momentum of the order of the saturation scale. This viewpoint is based on the fact that the system of partons 
(gluons) creates a new state of matter, the so-called Color Glass Condensate, in which the gluon density reaches the 
limited values of the order of l/a^ with new typical transverse momentum (the saturation scale). We developed a 
model that includes the gluon saturation and demonstrated that this model is able to describe both the inclusive 
hadron production at high energies including the first data from the LHC and the deep inelastic scattering data from 
HERA in a unique fashion. 

We predicted an increase of dNch/dy\ri=0: niean transverse momentum and the multiplicity of produced charged 
hadrons with energy which is in accordance with the first LHC data measured by ALICE [l|, |j] , CMS ^ and ATLAS 
collaboration, see Figs. I2I3I6I In the framework of high density QCD all these phenomena are closely related to 
the growth of the saturation momentum as a function of energy and of density of partons. It should be stressed 
that the other high-energy phenomenological approaches [5i| cannot describe the dependence of the average transverse 
momentum of the produced hadron on energy and hadron multiplicities. 

We showed that recently reported data by the CMS collaboration \2\ on the differential yield of charged hadrons 
at low pt for ^/s = 2.36 TeV reveal an interesting information on the mini-jets production and its connection with 
the saturation. We showed that the appearance of a peak in differential yield of charged hadrons at low pT is closely 
related to the mini-jet mass and the value of the saturation scale. 

We provided various predictions for the upcoming LHC measurements at higher energies in pp collisions. We believe 
that this paper will be useful for the microscopic interpretation of the upcoming LHC data and will lead to a deeper 
understanding of the hadron interactions at high energy in the framework of QCD. 

Concluding this paper we would like to answer the question: what can be here considered as a possible signal of the 
saturation (CGC) which are not contaminated with the non-perturbative physics related to unknown confinement of 
quarks and gluon? The main non-perturbative parameter that we have to introduce is rrijet ■ The rapidity distribution 
dNch/dr] at \r]\ < 1 (Fig. [5] ), the pt spectrum of hadron at low px < TTijet and the position of the maximum in 
d^Nch/drjdpT (Fig. [5]) depend on the value of rrijet and the success of our description indicates that we have chosen 
this parameter in self-consistent way. However, the energy dependence of dNch/drj at jr/j < 3.5 and the average value 
of the transverse momentum < px > oi hadrons as well as the multiplicity dependence of < pt > and the rapidity 
dependence of the maximum in d'^Nch/drjdpT for jryj < 3.5 are the typical consequences of the saturation approach 
since the main contribution in the calculations of these observables is originated from the transverse momenta of the 
order of Qs- Two factors determine the behavior of the observables at \r]\ > 3.5 : (1 — x)'^ suppression of the gluon 
densities in projectile and the increase of the saturation momentum in the target. Since the (1 — a;)^ factor reflects 
the well-known behavior of the structure function F2 at large-x, this factor will be the same in all other approaches 
while the additional increase due to the energy dependence is a typical features of the saturation approach. Notice 
also that at LHC energy ^/s — 14 TeV, the contribution of (1 — x)^ correction of unintegrated gluon density within 
the rapidities region considered here (Fig. [2]) is negligible and at 7 TeV this contribution is less than 5%. 

The above discussion shows that the comparison of our prediction with the high LHC energy data will be crucial 
for our approach. We are happy to make predictions before the experimental data from the LHC at high energy. We 
believe that if the coming data confirms our predictions, this will be indeed a first important step toward discovery 
of the CGC phase of the matter at LHC. The fact that we had to introduce several phenomenological parameters 
reflects our lack of theoretical knowledge for quark and gluon confinement and cannot be overcome in any models. 
Our experience tells us that when the data for higher energies will be published a lot of phenomenological models 
will appear but the CGC (saturation) approach is the only one that gives the predictions. It has happened once for 
nucleus-nucleus scattering at RHIC and, we hope that the situation will repeat itself at the LHC. 

The particle production scheme presented in this paper can be also applied to the calculation of inclusive hadron 
production in heavy ion collisions at LHC. We are currently working on this problem and plan to report on this in 
the near future. 
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